RECURSIVE DISCRETE FOURIER TRANSFORMATION APPARATUS 



BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an arithmetic operation 
apparatus for Fourier transformation and more particularly 
to an arithmetic operation apparatus for carrying out Fourier 
transformation with simple arithmetic processing so as to 
obtain an operation result of Fourier transformation in a short 
time. 

2. Description of the Related Art 

Conventionally, the Fourier transformation method used 
in frequency analysis and the like of time-series data stream 
has been employed in high compression encoding of audio signal 
and moving picture signal as well as spectrum analysis of 
signals in audio signal processing field, moving picture data 
processing field for medical system and the like. Further, 
the Fourier transformation method has been widely used as 
modulation technology and demodulation technology in 
communication field. 

According to the Fourier transformation method, data 
stream sampled as digital amount is treated as N groups (N 
is an integer, for example, 1024) and a time interval, in which 
the N-data stream exist, is assumed as a window period while 
that window period is assumed to be basic frequency. Then, 
signal components of data stream existing in the window period 
are obtained as real part and imaginary part of a harmonic 
signal of the basic frequency. 

Then , such data stream is handled as discrete data stream 
sampled at every predetermined period and the discrete data 
is subjected to Fourier transformation. The method is called 
discrete Fourier transformation (DFT) . The discrete Fourier 
transformation technology is employed as an analysis 



technology for control technology, in which, for example, a 
state of production process is obtained with the discrete data, 
the quality of the process is kept optimum by analyzing the 
obtained data and the rate of acceptable products in 
5 manufactured products is improved. Its application field has 
been widening every year . 

According to the discrete Fourier transformation 
technology, supplied signals are sampled at every 
predetermined time interval and a voltage value obtained by 

10 the sampling is obtained as sampled data. When the data stream, 
which is a group of the obtained data, is N data obtained in 

a predetermined time t , x(t), x(t+l), x(t+2), x(t+N-2), 

x(t+N-l), a value x (k, t) of the discrete Fourier 
transformation obtained with respect to the N data is defined 

15 in a following equation. 
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20 

As evident from this equation, the Fourier 
transformation is carried out on supplied data stream just 
such that its inherent basic function is convoluted in every 
point to be obtained and an operation for convoluting the basic 
25 function is conducted by a plurality of multiplication 
processings . 

When that multiplication processing is carried out with 
a special multiplying circuit or a digital signal processor 
(DSP), it has been known that a load upon hardware used for 
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arithmetic operation in the multiplying circuit or the DSP 
is increased tremendously. 

As for the load on this multiplication processing, a 
transformation equation expressed with the aforementioned 
5 equation (1) needs multiplications of 4 X n 2 times, so that 
in case where N is 1024, for example, the multiplication 
frequency is about 4,200,000 times. Therefore, the load of 
the arithmetic operation is very large and further, if the 
number N of points to be handled as the data stream is increased , 

10 the multiplication processing frequency is increased by second 
power, which is unfavorable in terms of actual field. 

Thus, as the discrete Fourier transformation equation, 
fast Fourier transformation (FFT) , whose operation efficiency 
is raised by deforming the stream by paying attention to the 

15 regularity that the basic function is composed of a cyclic 
function, is employed actually. 

An arithmetic operation method called butterfly 
operation is employed for the FFT. The butterfly operation 
is so constructed that a simple integer, for example, 2 is 

20 defined as a basic number and supplied binary complex data 
is subjected to complex operation comprised of a single 
addition, subtraction and multiplication each so as to output 
binary complex data. 

Therefore, the N-point FFT is comprised of log 2 N-step 

25 stage and (N/2)log 2 N butterfly operations and an operation 
result of the FFT can be obtained depending on the number of 
multiplications of (2N)log 2 N times. Consequently, the 
N-point FFT is used as a Fourier transformation method ensuring 
a high operating efficiency. 

30 Usually, the FFT or discrete Fourier transformation 

(DFT) is used in such a manner, so that data stream, sampled 
successively at a predetermined time interval and then supplied 
is subjected to Fourier transformation. In that Fourier 
transformation, the sampled data stream is stored in a memory 

35 circuit in succession temporarily and when the number of data 
stored in the memory circuit temporarily reaches N , the Fourier 
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transformation processing for that N data is started. 

In a period of the transformation processing , data stored 
temporarily in the memory circuit is kept as it was memorized 
and subjected to Fourier transformation processing. When the 
5 arithmetic operation is terminated, new data is supplied to 
the memory circuit and stored temporarily. After a 
predetermined amount of N is stored, next transformation 
processing is started. 

However, because data is supplied without any 

10 interruption even during Fourier arithmetic operation, 
according to another available method, a memory circuit for 
storing N data streams is provided as another system and each 
time when N data is supplied, temporary storage processing 
and arithmetic operation are carried out alternately and 

15 Fourier transformation processing is carried out continuously 
on time-series data. 

Because the Fourier transformation conducts arithmetic 
operation processing on N temporarily stored data streams while 
handling the data stream as a block unit, a delay time 

20 corresponding to at least N data is generated in such Fourier 
transformation processing, so that any transformation 
processing result cannot be obtained at real time. 

Thus, although data stream is supplied successively, 
a result of Fourier analysis cannot be obtained until each 

25 sampling of N data streams is finished. 

Contrary to this, in order to carry out Fourier 
transformation processing on N-point data streams containing 
new data supplied successively to be sampled and obtain a result 
of the Fourier transformation processing each time when sampled 

30 data is supplied, it comes that Fourier transformation 
processing of N points has to be carried out in a single sampling 
period and even in case where FFT arithmetic operation method 
developed for such high- speed arithmetic operation is employed, 
N-time arithmetic operation speed is demanded to obtain a 

35 result of Fourier transformation executed continuously from 
the FFT. However, in normal condition, it is difficult to 



execute such ultra FFT and supply a result of the arithmetic 
operation. 

In Fourier transformation, generally, FFT arithmetic 
operation processing method is employed. Usually, the FFT 
5 arithmetic operation processing is carried out on N-point data 
stream supplied by quant it izing sampling frequency with fs 
at a frequency interval of fs/N. 

Generally, even during Fourier transformation, data is 
fetched into other buffer memory so that the Fourier analysis 

10 processing is carried out continuously to time- series data 
supplied during arithmetic operation. In a period in which 
data temporarily stored in one buffer memory is being subjected 
to arithmetic operation processing, data supplied to the other 
buffer memory is fetched . When fetching of N data is terminated , 

15 the arithmetic operation and the data fetching processing may 
be switched to execute Fourier operation processing. However, 
in this case, two pairs of the buffer memories and FFT arithmetic 
operation processing means are required, which is not a 
favorable way in economic viewpoint. 

20 Further, because this method employs block processing 

method, in which N supplied time-series data are handled in 
batch, a result of the Fourier transformation to the fetched 
N-point data is not outputted until N sampling time passes. 
As an analysis result obtained at that time, only a result 

25 of the Fourier transformation processing on every N samples 
is outputted. 

That is, the result of the Fourier transformation on 
the newest N-point data containing data sampled successively 
and supplied cannot be outputted at real time. The Fourier 

30 transformation processing at every sampling period is 
necessary to output the result of the transformation at real 
time. However, if the Fourier transformation is carried out 
continuously at every sampling time interval, the amount of 
arithmetic operation per unit time is enormously increased, 

35 which is a not realistic method. 

On the other hand, as a method for obtaining a result 
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of continuous Fourier transformation processing, "Fourier 
transformation device and Fourier transformation method" has 
been proposed in Japanese Patent Application Laid-Open No. 
Hl-59454. 

5 This publication describes a method for 

Fourier-transforming a vibration waveform value sampled and 
supplied. According to the same method, a difference between 
a newly supplied vibration waveform value and an old vibration 
waveform value already supplied and used for Fourier 

10 transformation processing is obtained and a newer complex 
vibration value than the old complex vibration value provided 
by the arithmetic operation is obtained each time when a 
sampling value of the vibration waveform is supplied. 

Although as an application capable of carrying out the 

15 Fourier transformation continuously, frequency analysis may 
be considered, if such application is taken into account, any 
frequency band is required to be analyzed at any resolution. 

Japanese Patent Application Laid-Open No. H3-63875 has 
disclosed "A computation method for discrete Fourier 

20 transformation using cyclic technology" , in which discrete 
Fourier transformation for carrying out energy spectrum 
calculation with an ultrasonic diagnosis apparatus is 
described. The same publication describes discrete Fourier 
transformation, which continuously computes using a sample 

25 data file by cyclic filter operation. 

However, although this invention has described the 
method for achieving a circuit to be configured with complex 
data as input, it has not considered a cyclic Fourier 
transformation method having a simple structure , in which only 

30 real-part signals are supplied to transform signals supplied 
as time-series data for use in frequency analysis. 

SUMMARY OF THE INVENTION 

35 Accordingly, an object of the present invention is to 

provide a recursive discrete Fourier transformation device 
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in which Fourier analysis is configured for supplied real-part 
signals in a simple way, and more particularly, a recursive 
discrete Fourier transformation device suitable for frequency 
analysis to be executed at any desired resolution in any desired 
5 frequency band. 

To achieve the above object, there is provided a 
recursive discrete Fourier transformation device wherein data 
values x(t), x(t+l), x(t+2), x(t+3), . .., x(t+N-l), x(t+N) 
sampled at times t, t+l, t+2, t+3, t+N-1, t+N (N is a 

10 positive integer which is 1 or more) each having an equal 
interval are supplied and as complex Fourier coefficients under 
degree k (k is 0 or a positive integer smaller than N) obtained 
by, with such N data values supplied since time t as a data 
stream, carrying out complex Fourier transformation on the 

15 data stream, a real part X r (k, t) and an imaginary part Xi(k, 
t) are obtained, the discrete Fourier transformation device 
comprising: a first temporary storage means for storing the 
data stream x(t), x(t+l), x(t+2). x(t+3), x(t+N-l) 
supplied since time t at time t+N-1 temporarily; a discrete 

20 Fourier operation means for obtaining the complex Fourier 
coefficients X r (k, t) and Xi(k, t) of the data stream stored 
temporarily in the first storage means; and a second temporary 
storage means for storing the complex Fourier coefficients 
X r (k, t) andXi(k, t) obtained by the discrete Fourier operation 

25 means; the discrete Fourier operation means including: a 
subtracting portion for obtaining a data value of a difference 
between a data value x(t+N) supplied at time t+N and a data 
value x(t) memorized temporarily in the first storage means; 
a constant multiplying portion for obtaining a signal with 

30 a predetermined amplitude by multiplying the obtained data 
value of the difference with a positive constant value A for 
giving a predetermined amplitude; an adder portion for 
obtaining a summed signal by summing the signal with the 
predetermined amplitude obtained from the constant 

35 multiplying portion and one of the real part X r (k, t) and the 
imaginary part X±(k, t) of the complex Fourier coefficients 
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stored temporarily in the second temporary storage means ; and 
a basic function arithmetic processing portion for receiving 
the summed signal obtained from the adder portion and the other 
of the real part X r (k, t) and the imaginary part X±(k, t) of 
5 the complex Fourier coefficients stored temporarily by the 
second temporary storage means and carrying out an arithmetic 
operation on the received signals using a constant based on 
a basic frequency thereby to obtain the complex Fourier 
coefficients X r (k, t+1) and Xj.(k, t+1) at time t+1. 

10 According to the present invention, data stream 

comprised of supplied N samples is subjected to complex Fourier 
transformation and then, the oldest sample value used for 
complex Fourier transformation is subtracted from the value 
of a new sample to be supplied next so as to obtain a subtraction 

15 value . By deleting the oldest sample at the same time , discrete 
Fourier arithmetic operation result on new N samples, which 
are renewed by each sample, can be obtained based on that 
subtraction value and the result of the complex Fourier 
operation. Thus, different from the conventional Fourier 

20 operation which is carried out after N sample data is supplied, 
it is possible to configure a recursive Fourier transformation 
device, which obtains a result of fast Fourier operation on 
points N so as to provide with an analysis resolution necessary 
in a sample period. 

25 According to a preferred embodiment of the present 

invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) can be set selectively with 1, square root of 
N or l/N. 

30 According to this embodiment, in addition to the above 

described effect, the operation can be carried out by FFT in 
which the positive constant value A is set up by selecting 
appropriately 1, N or a square root of N. This is employed 
in, for example, communication apparatus. In case where 

35 information signal to be transmitted is subjected to inverse 
fast Fourier transformation (IFFT) processing to be converted 



to time-series signal and that time-series signal is 
transmitted, according to Fourier operation method for 
obtaining information signal by supplying that transmitted 
signal to the FFT, which executes the FFT operation, the 
5 constant value A corresponding to the constant value a, which 
is, for example, 1, N or the square root of N used for the 
IFFT operation, is employed in the FFT operation processing 
circuit, so as to construct the FFT which operates in a 
complementary way with the IFFT to decode the information 

10 signal. By using the constant value A corresponding to the 
value a, used for the IFFT operation for a signal supplied 
from other mating system, a recursive Fourier transformation 
device having an excellent characteristic can be configured. 
Further, to achieve the above object, there is provided 

15 a recursive discrete Fourier transformation device wherein 
datavaluesx(t) , x(t+l) , x(t+2) , x(t+3) , . . . , x(t+N-l) , x(t+N) 
sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 
interval are supplied and as complex Fourier coefficients under 

20 degree k (k is 0 or a positive integer smaller than N) obtained 
by, with such N data values supplied since time t as a data 
stream, carrying out complex Fourier transformation on the 
data stream, a real part X r (k, t) and an imaginary part X ± (k, 
t) are obtained, the discrete Fourier transformation device 

25 comprising: a first temporary storage means for storing the 
data stream x(t), x(t+l), x(t+2), x(t+3), x(t+N-l) 
supplied since time t at time t+N-1 temporarily; a discrete 
Fourier operation means for obtaining the complex Fourier 
coefficients X r (k, t) and X ± (k, t) of the data stream stored 

30 temporarily in the first storage means ; and a second temporary 
storage means for storing the complex Fourier coefficients 
X r (k, t) andX ± (k, t) obtained by the discrete Fourier operation 
means , wherein the discrete Fourier operation means obtains 
complex Fourier coefficients X r (k, t) and Xi(k. t) according 

35 to following equations. 
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X, (k,t + 1) - \x r (k,l) + ±[x(t + JV)- *(()]}xcos|2^] + *,(*,, )rin 

X,(t,, +1 )-Jr,.(t, ( )cos[2^]-j^(v) + iW' + W)-^)]j si „ [ 2 £] 

Where, A is a positive constant value for providing 
[x( t+N) -x(t ) ] with an amplitude. 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
transformation and then, by carrying out arithmetic operation 
according to a predetermined equation based on a value of a 
new sample supplied next, the oldest sample value used for 
complex Fourier transformation and a result of the already 
obtained complex Fourier operation result, discrete Fourier 
operation result on new N samples, which are renewed by each 
sample, can be obtained. Thus, different from the 
conventional Fourier operation which is carried out after N 
sample data are supplied, it is possible to configure a 
recursive Fourier transformation device, which obtains a 
result of fast Fourier operation on points N so as to provide 
with an analysis resolution necessary in a sample period. 

According to a preferred embodiment of the present 
invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x( t ) is capable of being set selectively with 1 , square 
root of N or 1/N. 

Further, to achieve the above object, there is provided 
a recursive discrete Fourier transformation device wherein 
datavaluesx(t) , x(t + l) , x(t+2) , x(t+3) , . ... , x(t+N-l) , x(t+N) 

sampled at times t, t+1, t+2, t+3, , t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
since time t as a data stream, complex Fourier transformation 
is carried out to the data stream using a plurality of degrees 
k (k is 0 or a positive integer smaller than N) so as to obtain 
real parts X r (k, t) and imaginary parts Xi(k, t) as plural 
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sets of complex Fourier coefficients, the discrete Fourier 
transformation device comprising: a first temporary storage 
means for storing the data stream x(t), x(t+l), x(t+2), 
x(t+3), x(t+N-l) supplied since time t at time t+N-1 

5 temporarily; plural discrete Fourier operation means for 
obtaining the complex Fourier coefficients X r (k, t) and X ± (k, 
t) for the data stream stored temporarily in the first storage 
means for each of plural k values ; and a second temporary storage 
means for storing each set of the complex Fourier coefficients 

10 X r (k, t) and Xi(k, t) obtained by the plural discrete Fourier 
operation means corresponding to each k value, the discrete 
Fourier operation means including: a subtracting portion for 
obtaining a data value of a difference between a data value 
x(t+N) supplied at time t+N and a data value x(t) memorized 

15 temporarily in the first storage means ; a constant multiplying 
portion for obtaining a signal with a predetermined amplitude 
by multiplying the data value of the difference obtained by 
the subtracting portion with a positive constant value A for 
giving a predetermined amplitude; an adder portion for 

20 obtaining a summed signal by summing the signal with the 
predetermined amplitude obtained from the constant 
multiplying portion and one of a real part X r (k, t) and an 
imaginary part (k, t) of the complex Fourier coefficients 
stored temporarily by the second temporary storage means; and 

25 a basic function arithmetic processing portion for receiving 
the summed signal obtained from the adder portion and the other 
of the real part X r (k, t) and the imaginary part X±(k, t) of 
the complex Fourier coefficients stored temporarily in the 
second temporary storage means and carrying out an arithmetic 

30 operation on the received signals using a constant based on 
a basic frequency thereby to obtain complex Fourier 
coefficients X r (k, t+1) and X ± (k, t+1) at time t + 1. 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 

35 transformation in terms of plural basic frequencies and then, 
the oldest sample value used for complex Fourier transformation 
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is subtracted from the value of a new sample to be supplied 
next so as to obtain a subtraction value. By deleting the 
oldest sample at the same time, discrete Fourier operation 
result on new N samples, which are renewed by each sample, 
5 can be obtained based on that subtraction value and the result 
of the complex Fourier operation in terms of plural basic 
frequencies. Thus, different from the conventional Fourier 
operation which is carried out after N sample data is supplied, 
it is possible to configure a recursive Fourier transformation 
10 device, which is operated fast to obtain a result of Fourier 
operation in terms of plural basic frequencies within a sample 
period. 

According to a preferred embodiment of the present 
invention, the quantity of the degrees k is N. 
15 According to this embodiment, in addition to the 

above -described effect, a recursive Fourier transformation 
device, which is operated fast so as to obtain a result of 
Fourier operation for all N basic frequencies, can be 
configured. 

20 According to a preferred embodiment of the present 

invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) is capable of being set selectively with 1, square 
root of N or 1/N. 

25 Further, to achieve the above object, there is provided 

a recursive discrete Fourier transformation device wherein 

datavaluesx(t) , x(t+l) , x(t+2) , x{t+3) , x(t+N-l) , x(t+N) 

sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 

30 interval are supplied and with such N data values supplied 
since time t as a data stream, complex Fourier transformation 
is carried out to the data stream using a plurality of degrees 
k (k is 0 or a positive integer smaller than N) so as to obtain 
real parts X r (k, t) and imaginary parts Xi(k, t) as plural 

35 sets of complex Fourier coefficients, the discrete Fourier 
transformation device comprising: a first temporary storage 
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means for storing the data stream x(t), x(t+l), x(t+2), 
x(t+3), x(t+N-l) supplied since time t at time t+N-1 

temporarily; plural discrete Fourier operation means for 
obtaining the complex Fourier coefficients X r (k, t) and Xi(k, 
5 t ) for the data stream stored temporarily in the first storage 
means for each k value; and a second temporary storage means 
for storing each set of complex Fourier coefficients X r (k, 
t) and Xi(k, t) obtained by the plural discrete Fourier 
operation means corresponding to each k value, the discrete 

10 Fourier operation means including: a common subtracting 
portion for obtaining a data value of a difference between 
a data value x(t+N) supplied at time t+N and a data value x(t) 
memorized temporarily in the first storage means; a common 
constant multiplying portion for obtaining a signal with a 

15 predetermined amplitude by multiplying the data value of the 
difference obtained by the common subtracting portion with 
a positive constant value A for giving a predetermined 
amplitude; an adder portion for obtaining a summed signal by 
summing the signal with the predetermined amplitude obtained 

20 from the common constant multiplying portion and one of a real 
partX r (k, t) and an imaginary part (k, t) of the complex Fourier 
coefficients stored temporarily in the second temporary 
storage means; and a basic function arithmetic processing 
portion for receiving the summed signal obtained from the adder 

25 portion and the other of the real part X r (k, t ) and the imaginary 
part Xi(k, t) of the complex Fourier coefficients stored 
temporarily in the second temporary storage means and carrying 
out an arithmetic operation on the received signals using a 
constant based on a basic frequency thereby to obtain the 

30 complex Fourier coefficients X r (k, t+1) andX ± (k, t+1) at time 
t+1. 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
transformation in terms of plural basic frequencies and then, 
35 the oldest sample value used for complex Fourier transformation 
is subtracted from the value of a new sample to be supplied 
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next so as to obtain a subtraction value. By deleting the 
oldest sample at the same time, discrete Fourier operation 
result on new N samples, which are renewed by each sample, 
can be obtained based on that subtraction value and the result 
5 of the complex Fourier operation in terms of plural basic 
frequencies. At this time, the first primary storage means, 
subtracting portion and constant multiplying portion can be 
used in common . Thus , different from the conventional Fourier 
operation which is carried out after N sample data is supplied, 

10 a recursive Fourier transformation device, which is operated 
fast to obtain a result of Fourier operation in terms of plural 
basic frequencies within a sample period, can be configured. 

According to a preferred embodiment of the present 
invention, the quantity of the degrees k is N. 

15 According to a preferred embodiment of the present 

invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) is capable of being set selectively with 1, square 
root of N or 1/N. 

20 Further, to achieve the above object, there is provided 

a recursive discrete Fourier transformation device wherein 

datavaluesx(t) , x(t + l) , x(t+2) , x(t + 3) , , x(t+N-l) , x(t+N) 

sampled at times t, t+1, t+2, t+3, .... t+N-1, t+N (N is a 
positive integer which is 1 or more) each having an equal 

25 interval are supplied and as a complex Fourier coefficient 
under degree k (k is 0 or a positive integer smaller than N) 
obtained by, with such N data values supplied since time t 
as a data stream, carrying out complex Fourier transformation 
on the data stream, a real part X r (k, t) and an imaginary part 

30 X ± (k, t) are obtained, the discrete Fourier transformation 
device comprising: a data updating means for obtaining a first 
subtraction signal by subtracting data x(t) supplied before 
N sampling period from data x(t+N) supplied at time t+N; a 
recursive processing means for obtaining a new second 

35 subtraction signal by subtracting an addition signal generated 
recursively using an already generated second subtraction 
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signal from the obtained first subtraction signal; and a 
multiplying means for obtaining the real part X r (k, t) of the 
Fourier coefficients by summing up a signal obtained by 
multiplying the new second subtraction signal obtained by the 
5 recursive processing means with a first constant value and 
a signal obtained by multiplying the second subtraction signal 
supplied before a sampling period with a second constant value 
and for obtaining the imaginary part Xi(k, t) of the Fourier 
coefficients by multiplying the new second subtraction signal 

10 with a third constant value, wherein the addition signal 
generated recursively by the recursive processing means is 
a signal obtained by summing up a signal obtained by multiplying 
the second subtraction signal obtained before a sampling period 
with a fourth constant value and the second subtraction signal 

15 obtained before two sampling periods . 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
transformation and then, the oldest sample value used for 
complex Fourier transformation is subtracted from the value 

20 of a new sample to be supplied next so as to obtain a subtraction 
value . By deleting the oldest sample at the same time , discrete 
Fourier operation result on new N samples, which are renewed 
by each sample , can be obtained based on that subtraction value 
and the result of the complex Fourier operation. Thus, 

25 different from the conventional Fourier operation which is 
carried out after N sample data is supplied , a recursive Fourier 
transformation device, which obtains a result of fast Fourier 
arithmetic operation on points N so as to provide with an 
analysis resolution necessary in a sample period can be 

30 configured. 

Further, to achieve the above object, there is provided 
a recursive discrete Fourier transformation device wherein 
datavaluesx(t) , x(t+l) , x(t+2) , x(t+3) , . . . , x(t+N-l) , x(t+N) 
sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N is a 

35 positive integer which is 1 or more) each having an equal 
interval are supplied and as a complex Fourier coefficient 
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under degree k (k is 0 or a positive integer smaller than N) 
obtained by, with such N data values supplied since time t 
as a data stream, carrying out complex Fourier transformation 
on the data stream, a real part X r (k, t) and an imaginary part 
5 X±(k, t) are obtained, the discrete Fourier transformation 
device comprising: a data updating means for obtaining a first 
subtraction signal by subtracting data x(t) supplied before 
N sampling period from data x(t+N) supplied at time t+N; a 
recursive processing means for obtaining a new second 

10 subtraction signal by subtracting an addition signal generated 
recursively using an already generated second subtraction 
signal from the obtained first subtraction signal ; and a 
multiplying means for obtaining the real part X r (k, t) of the 
Fourier coefficients by summing up a signal obtained by 

15 multiplying the new second subtraction signal obtained by the 
recursive processing means with a first constant value and 
a signal obtained by multiplying the second subtraction signal 
supplied before a sampling period with the second constant 
and for obtaining the imaginary part Xi(k, t) of the Fourier 

20 coefficients by multiplying the new second subtraction signal 
with a third constant value, wherein a transfer function H(Z) 
for the data updating means, the recursive processing means 
and the multiplying means connected as subsidiary components 
is given according to a following equation. 



cos 2 — - /sin 2 
1 N \ J I N 



where A is a positive constant value for providing 
[x(t+N)-x(t) ] with an amplitude. 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
30 transformation and then, by carrying out arithmetic operation 
according to a predetermined equation based on a value of a 
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new sample supplied next, the oldest sample value used for 
complex Fourier transformation and a result of the already 
obtained complex Fourier operation result, discrete Fourier 
operation result on new N samples, which are renewed by each 
5 sample , can be obtained . Thus , different from the conventional 
Fourier operation which is carried out after N sample data 
is supplied, a recursive Fourier transformation device, which 
obtains a result of fast Fourier arithmetic operation on points 
N so as to provide with an analysis resolution necessary in 

10 a sample period can be configured. 

According to a preferred embodiment of the present 
invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) is capable of being set selectively with 1, an 

15 inverse number of square root of N or 1/N. 

According to this embodiment, in addition to the above 
described effect, the operation can be carried out by FFT in 
which the positive constant value A is set up by selecting 
appropriately 1, N or an inverse number of square root of N. 

20 This is employed in, for example, communication apparatus. 
In case where information signal to be transmitted is subjected 
to inverse fast Fourier transformation (IFFT) processing to 
be converted to time-series signal and that time-series signal 
is transmitted, in Fourier operation method for obtaining 

25 information signal by supplying that transmitted signal to 
the FFT for executing the FFT operation, the constant value 
A corresponding to the constant value a, which is , for example , 
1, N or an inverse number of the square root of N used for 
the IFFT operation, is employed in the FFT operation processing 

30 circuit , so as to construct the FFT which operates in a 
complementary way with the IFFT to decode the information 
signal. By using the constant value A corresponding to the 
value a used for the IFFT operation for a signal supplied from 
other mating system, a recursive Fourier transformation device 

35 having an excellent characteristic can be configured. 

Further, to achieve the above object, there is provided 
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a recursive discrete Fourier transformation device wherein 
datavaluesx(t) , x(t + l) , x(t+2) , x(t + 3) , . . . , x(t+N-l) , x(t+N) 
sampled at times t, t + 1, t+2, t + 3, .... t+N-1, t+N (N is a 
positive integer which is 1 or more) each having an equal 
5 interval are supplied and with such N data values supplied 
since time t as a data stream, complex Fourier transformation 
is carried out to the data stream using a plurality of degrees 
k (k is 0 or a positive integer smaller than N) so as to obtain 
real parts X r (k, t) and imaginary parts X ± (k, t) as plural 

10 sets of complex Fourier coefficients, the discrete Fourier 
transformation device comprising: plural data updating means 
corresponding to the plurality of degrees k, for obtaining 
a first subtraction signal by subtracting data x(t) supplied 
before N sampling period from data x(t+N) supplied at time 

15 t+N; plural recursive processing means corresponding to the 
plurality of degrees k, for obtaining a new second subtraction 
signal by subtracting an addition signal generated recursively 
using an already generated second subtraction signal from the 
obtained first subtraction signal; and plural multiplying 

20 means corresponding to the plurality of degrees k, for 
obtaining a real part X r (k, t) of the Fourier coefficients 
by summing up a signal obtained by multiplying the new second 
subtraction signal obtained by the recursive processing means 
with a first constant value and a signal obtained by multiplying 

25 the second subtraction signal supplied before a sampling period 
with the second constant and for obtaining an imaginary part 
Xi(k, t) of the Fourier coefficients by multiplying the new 
second subtraction signal with a third constant value , wherein 
the addition signal generated recursively by each of the plural 

30 recursive processing means is a signal obtained by summing 
up a signal obtained by multiplying the second subtraction 
signal obtained before a sampling period with a fourth constant 
value corresponding to each degree k, and the second 
subtraction signal obtained before two sampling periods. 

35 According to the present invention, data stream 

comprised of supplied N samples is subjected to complex Fourier 
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transformation in terms of plural basic frequencies and then, 
the oldest sample value used for complex Fourier transformation 
is subtracted from the value of a new sample to be supplied 
next so as to obtain a subtraction value. By deleting the 
5 oldest sample at the same time, discrete Fourier operation 
result on new N samples, which are renewed by each sample, 
can be obtained based on that subtraction value and the result 
of the complex Fourier operation in terms of plural basic 
frequencies. Thus, different from the conventional Fourier 

10 operation which is carried out after N sample data is supplied, 
a recursive Fourier transformation device, which is operated 
fast to obtain a result of Fourier operation in terms of plural 
basic frequencies within a sample period, can be configured. 

According to a preferred embodiment of the present 

15 invention, the quantity of the degrees k is N. 

According to this embodiment, in addition to the above 
described effect, a recursive Fourier transformation device, 
which is operated fast so as to obtain a result of Fourier 
operation for all N basic frequencies, can be configured. 

20 Further, to achieve the above object, there is provided 

a recursive discrete Fourier transformation device wherein 

datavaluesx(t) , x(t+l) , x{t+2) , x(t+3) , x(t+N-l) , x(t+N) 

sampled at times t, t+1, t+2, t+3, , t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 

25 interval are supplied, data x(t) supplied before N sampling 
period is subtracted from data x(t+N) supplied at time t+N 
so as to obtain a first subtraction signal, and with such N 
data values supplied since time t as a data stream based on 
the obtained first subtraction signal, a complex Fourier 

30 transformation is carried out to the data stream using a 
plurality of degrees k (k is 0 or a positive integer smaller 
than N) so as to obtain real parts X r (k, t) and imaginary parts 
X ± (k, t) as plural sets of complex Fourier coefficients , the 
discrete Fourier transformation device comprising: plural 

35 recursive processing means corresponding to the plurality of 
degrees k, for obtaining a new second subtraction signal by 
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subtracting an addition signal generated recursively using 
an already generated second subtraction signal from the 
obtained first subtraction signal; and plural multiplying 
means corresponding to the plurality of degrees k, for 
5 obtaining a real part X r (k, t) of the Fourier coefficients 
by summing up a signal obtained by multiplying the new second 
subtraction signal obtained by the recursive processing means 
with a first constant value and a signal obtained by multiplying 
the second subtraction signal supplied before a sampling period 

10 with the second constant and for obtaining an imaginary part 
X ± (k, t) of the Fourier coefficients by multiplying the new 
second subtraction signal with a third constant value , wherein 
the addition signal generated recursively by each of the plural 
recursive processing means is a signal obtained by summing 

15 up a signal obtained by multiplying the second subtraction 
signal obtained before a sampling period with a fourth constant 
value corresponding to each degree k, and the second 
subtraction signal obtained before two sampling periods. 

According to a preferred embodiment of the present 

20 invention, the quantity of the degrees k is N. 

The nature, principle and utility of the invention will 
become more apparent from the following detailed description 
when read in conjunction with the accompanying drawings. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

In the accompanying drawings: 

Fig.l is a diagram showing schematically a recursive 
discrete Fourier transformation device according to a first 
30 embodiment of the present invention; 

Fig. 2 is a diagram showing schematically a relation 
between data value obtained by sampling the waveform of signal 
supplied to the recursive discrete Fourier transformation 
device of the first embodiment of the present invention and 
35 a corresponding DFT operation; 

Fig . 3 is a diagram showing the structure of the recursive 
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discrete Fourier transformation device according to the first 
embodiment of the present invention; 

Fig. 4 is a diagram showing in detail the structure of 
the recursive discrete Fourier transformation device 
5 according to the first embodiment of the present invention; 

Fig . 5 is a diagram showing the structure of the recursive 
discrete Fourier transformation device which computes Fourier 
coefficient with respect to basic frequency of N points 
according to the first embodiment of the present invention; 
10 Fig • 6 is a diagram showing the structure of the recursive 

discrete Fourier transformation device which computes Fourier 
coefficient with respect to basic frequency of N points 
according to the first embodiment of the present invention; 

Fig. 7 is a diagram showing a schematic structure of the 
15 recursive discrete Fourier transformation device according 
to a second embodiment of the present invention; 

Fig . 8 is a diagram showing the structure of the recursive 
discrete Fourier transformation device which carries out 
arithmetic operation using a transfer function according to 
20 the second embodiment of the present invention; 

Fig. 9 is a diagram showing other structure of the 
recursive discrete Fourier transformation device which 
carries out arithmetic operation using a transfer function 
according to the second embodiment of the present invention; 
25 and 

Fig. 10 is a diagram showing the structure of the 
recursive discrete Fourier transformation device which 
computes Fourier coefficient with respect to the basic 
frequency of N points according to the second embodiment of 
30 the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Hereinafter, the preferred embodiment of the recursive 
35 discrete Fourier transformation device will be described with 
reference to the accompanying drawings. 
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< First Embodiment > 

Fig . 1 shows schematically the structure of the recursive 
discrete Fourier transformation device according to the first 
embodiment of the present invention, which will be described 
5 with the same Figure. 

The recursive discrete Fourier transformation device 
comprises a data updating portion 1 to which data sampled in 
a predetermined time interval is supplied and which stores 
temporarily supplied new N (N: positive integer) data therein, 
10 a basic frequency settingportion 2 in which the basic frequency 
for carrying out the discrete Fourier transformation is set 
up, a recursive DFT operating portion 3 for carrying out 
recursive DFT operation, and a memory 4 for storing temporarily 
computed data. 

15 Next, an operation of the recursive discrete Fourier 

transformation device having such a structure will be 
described. 

First, supplied data is sampled with a sampling circuit 
(not shown) every predetermined time interval and then, the 

20 sampled and quantitized discrete data is supplied to the data 
updating portion 1 . 

The sampling circuit samples data values supplied at 
every time of a predetermined interval, t, t+1, t+2, t + 3, 
t+N-1, t+N (N: natural number ) and generates data streams x(t) , 

25 x(t+l), x(t+2), x(t+3), x(t+N-l), x(t+N) by converting 

the data values supplied at the respective times to sampling 
values corresponding to the respective times . 

The operation of the sampling circuit is the same as 
that of an A/D converter which converts supplied analog signal 

30 to digital signal, converting an analog signal voltage value 
supplied at a predetermined time interval, given in the form 
of an inverse number of a sampling frequency, to a digital 
signal value . The converted digital signal value is a voltage 
having an analogous relation with an analog voltage, which 

35 gives a pulse amplitude modulation signal or that voltage value 
is expressed with a binary digital value. 
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In this manner, sampling data x(t) sampled at time t, 
x(t+l) sampled at time t + l # . . . , -x(t+N-l) sampled at t+N-1, 
x(t+N) sampled at time t+N, .... are supplied to the data 
updating portion 1. 
5 The data-updating portion 1 updates the newest N data 

(N: positive integer) of the supplied data, and stores them 
temporarily. 

That is, data supply is started from x(t) and when data 
x(t+l), x(t+2), x(t+3) are supplied, all the supplied data 

10 x(t), x(t+l), x(t+2), x(t+3) are stored temporarily. This 
temporary storage operation is continued until data x( t+N-1) 
is entered and when data x{ t+N-1) is entered, the total data 
number becomes N, so that data region of the data updating 
portion becomes full. 

15 When next data x( t+N) is supplied with such a condition, 

the total data number becomes N+l. Thus, the data updating 
portion 1 subtracts x(t) from x(t+N) and then supplies data 
obtained by that subtraction to the recursive DFT operating 
portion 3 while removing the oldest data x(t) from the memory 

20 portion. Consequently, the data updating portion 1 
temporarily stores Ndatax(t+1) , x(t+2) , x(t + 3) , . . . , x(t+N-l) 
x(t+N) . 

Likewise, when next datax(t+N+l) is supplied, the data 
updating portion 1 obtains x(t+N+l)-x(t+l) and supplies it 

25 to the recursive DFT operating portion 3 while removing x(t + l) 
from the memory. Consequently, the newest N data among 
supplied data is always stored in the data-updating portion 
1 temporarily. 

Such temporarily stored N data is supplied to the 

30 recursive DFT operating portion 3 and a last FFT operation 
result stored temporarily in the memory 4 is supplied to the 
recursive DFT operating portion 3 as recursive data. Then, 
the recursive DFT operating portion 3 carries out recursive 
discrete Fourier transformation based on a method described 

35 later, according to frequency resolution information set up 
by the basic frequency setting portion 2 and outputs a result 
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of the operation. 

Next, this recursive discrete Fourier transformation 
operation method will de described in detail with reference 
to a conventional Fourier operation as well. 
5 Fig. 2 explains schematically a relation between data 

value obtained by sampling the waveform of supplied signal 
in a sampling period ts and corresponding DFT operation. 

The same Figure shows a set of N real number data values 

x(t), x(t+l), x(t+2), x(t+3) x( t+N- 1) sampled in a 

10 predetermined sampling period from time t and a set of N real 
number data values x(t+l), x(t+2), x(t+3), x(t+N-l), 
x(t+N) sampled from time t+1. 

Then, the value X(k, t) of discrete Fourier 
transformation obtained with respect to data stream which is 
15 N real number data values sampled since time t is defined in 
the following expression. 

^(*.')-^"|^(«)exp|-/2(„-<)^j - (2) 

^P WroS [ 2( "-' ) f]-^'f Wsin [ 2( "- !) f] 

= X r (k,t)- jxXk,t) 
20 for k = 0, 1, , N-1 

Transformation of real part X r (k, t) and imaginary part 
X ± (k, t) of complex Fourier transformation of sampling data 
stream x(t) ~ x(t+N-l) supplied since arbitrary time t is 
25 defined as follows. 

Xrijl,t) = lW g^") cos [ 2 ("-0^] - O) 

xXkj) = ^p(n)sm^2(n-t)fj . - (4) 
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for k = 0, 1, ...,N-1 



The transformation equation for the sampling data 
supplied since time t is defined in the above manner. Next, 
5 a transformation equation for a sampling data stream to be 
supplied since time t+1 will be described below. 

That is, the sampling data stream supplied since the 
time t + 1 is expressed as a set of x(t+l), x(t+2) , x(t+3) , 
•x(t+N-l), x(t+N). Because data x(t) is deleted from data 
10 stream supplied since time t while new data is added, 
transformation of the real part is expressed in a following 
equation ( 5 ) and that expressed equation ( 5 ) is developed as 
follows. 




^(M + l) = ^|^)su,[2(„-,)^] - (7) 

20 Then, if x(t), which is the oldest data is deleted and 

x(t+N) which is the newest data is incorporated, a Fourier 
coefficient Y r (k, t+1) can be expressed with Fourier 
coefficient X r (k, t) and likewise, Yi(k, t+1) can be expressed 
with Xi(k, t) . 

25 As a result, the equations (6), (7) are expressed as 

follows . 
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Y r (k,t+l)-X r (k,t) + -j=[x(t + N)-x{t)] - (8) 

Y i (k t t + \)-'X i (k,t) - (9) 

Next , the transformation equation for the imaginary part 
is also developed in the same way. 
5 That is, when new data x(t+N) is supplied, the 

transformation of the imaginary part is also developed as shown 
in a following equation. 

Ar,(jk,*+l) = --^r^jc(/i)siii|2(/i-f-l)^-J — (10) 

10 = Y;.(*,f +l)cos -Y r (k,t +l)sin |2^j 

Therefore, the equations (5) , (10) can be expressed based 
on the relation of the equations (8), (9). 

X. t (k,t + 1) = jjT r (k,t) + -j=[x(t + N) - x(t)]x cos + X, (Msin|2^j|(ll) 

X, (k,t + 1) = X t {k,t)cos - |z r (k,t) + -j^Mt + N) - *(* )]j sin ( 12 ) 

15 In this way, the discrete Fourier transformation at time 

t + 1 can be executed according to the equations (11) , (12) by 
using the result of the discrete Fourier arithmetic operation 
at time t recursively and further using data value which is 
a difference between newly fetched sample value x(t+N) and 

20 a sample value x(t) to be deleted. 

Although the discrete Fourier transformation obtained 
here is FIR filter (acyclic digital filter) which originally 
carries out processing for finite sample number N , the discrete 



Fourier transformation at current sampling time is realized 
using IIR filter (cyclic digital filter), which introduces 
the discrete Fourier transformation using a result of the 
discrete Fourier transformation obtained at a previous 
5 sampling time, as evident from the aforementioned equations 
(11), (12). 

According to the method for achieving this DFT with the 
IIR filter, the hardware configuration for arithmetic 
operation can be made simpler than constructing the DFT with 

10 the FIR filter. Thus, if the recursive discrete Fourier 
transformation expressed with the aforementioned 
transformation equation is carried out, it is possible to 
supply new sampled data successively at every sampling time 
interval shown in Fig. 2 and carry out high efficiency Fourier 

15 transformation to the newest data stream containing new data 
values . 

In this way, the Fourier operation is enabled using the 
basic frequency expressed with cos(2jtk/N) and sin cos(2rck/N). 
If the basic frequency is assumed to be Tc, Ts and l/( square 
20 root of N) indicating the amplitude is assumed to be A, the 
equations (11), (12) are expressed as follows. 

X r (k,t + l)={X r (k,t) + A[x{t + N)-x{t)}Tc + X l {k,t)Ts - (13) 
X i (A:,? + l) = Z l .(^)rc-{Z r (^) + 4^ + ^")-^W}r5 - (14) 

Where , Tc = cos \l— 1 . Ts = sin \l— 1 , A=-^= 

25 In this way, the Fourier transformation operation is 

enabled using the equations (13), (14). Next, the structure 
of the recursive discrete Fourier transformation device , which 
achieves this, will be described. 

That is, the equation (13) indicates that the real part 

30 X r (k, t+1) of Fourier coefficient at time t+1 can be obtained 
at object basic frequencies To, Ts by using data x(k, t+N) , 
x(k, 1) supplied at time t+N and time t and real part X r (k, 
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t) and imaginary part Xi(k, t) of Fourier coefficient obtained 
at time t recursively. 

Then, the equation (14) indicates that the imaginary 
part Xi(k, t+1) of Fourier coefficient at time t+1 can be 
5 obtained. Next , a recursive discrete Fourier transformation 
device, which operates corresponding to these arithmetic 
operation equations, will be described. 

Fig. 3 is a recursive discrete Fourier transformation 
circuit composing major components of the recursive discrete 

10 Fourier transformation device and Fig. 4 is a diagram showing 
the basic function arithmetic processing portion 30 further 
in detail. The recursive discrete Fourier transformation 
device will be explained with reference to these Figures . 

Referring to Fig. 4, the recursive discrete Fourier 

15 transformation circuit comprises a delay circuit 11 and a 
subtractor 12, which compose a data updating portion 1, a 
plurality of multipliers 31-35, which compose a recursive DFT 
operating portion 3, a plurality of adders 36-38 and delay 
circuits 41, 42 which compose a memory portion 4. 

20 Next, an operation of the recursive discrete Fourier 

transformation circuit having such a structure will be 
described. 

First, supplied input signals are sampled at times x(t) , 
x(t + l) , x(t + 2) , x(t + 3) , . . . , x(t+N-l) , x(t+N) (N is a natural 
25 number ) with an equal interval by a sampling circuit ( not shown ) 
and then, discrete data that input signal levels at those times 
are gathered in the form of data stream x(t) , x(t+l) , x(t+2) , 
x(t+3) , . . . , x(t+N-l) , x(t+N) as sampling values are supplied 
to the data updating portion 1. 
30 In the data updating portion 1, one part of supplied 

data is supplied to a delay unit 11 while other part is supplied 
to the subtractor 12. The delay unit 11 delays a signal by 
a time corresponding to the period of N sampling clocks. 

That is, the delay unit 11 has a function for updating 
35 N data supplied newly while storing them temporarily. 

Here, data supply is started from x(t) and when data 
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x(t+l) , x(t+2) , x(t+3) are supplied, the delay unit 11 stores 
allthesupplieddatax(t) , x(t+l) , x(t+2) , x(t+3) temporarily. 
This temporary storage operation is continued until data 
x(t+N-l) is entered and when data x(t+N-l) is entered, the 
5 total data number becomes N, so that data region of the data 
updating portion becomes full. 

Because, when next data x(t+N) is supplied under this 
condition, total data quantity is N+l, the delay unit 11 
supplies the oldest data x(t) to the subtractor 12 and the 

10 subtractor 12 subtracts that oldest data x(t) from supplied 
datax(t+N). Then, [x(t+N) -x(t) ] obtained by the subtraction 
is supplied to the recursive DFT operating portion 3. 

The recursive DFT operating portion 3 carries out 
multiplication and addition of such supplied signals. 

15 First , the multiplier 31 multiplies the supplied signal 

[x( t+N) -x(t ) ] with the aforementioned constant A and the 
multiplied signal is supplied to the adder 36. 

A signal supplied from the multiplier 31 and a signal 
obtained by delaying the real part of an output signal from 

20 the recursive discrete Fourier operation circuit by a sampling 
period by the delay unit 41 are supplied to the adder 36 so 
that they are summed up and that summed signal is supplied 
to the multipliers 32, 33. 

In one multiplier 32, the supplied signal is multiplied 

25 with the basic frequency Tc = cos ( 2jtk/N) and then, a signal 
obtained by the multiplication is supplied as an input signal 
to the adder 37 . In another multiplier 33 , the supplied signal 
is multiplied with the basic frequency -rs = -sin( 2rck/N) and 
a signal obtained by that multiplication is supplied as one 

30 input signal to the adder 38. 

In other input signal to be supplied to the 
aforementioned adder 37, the imaginary part of an output signal 
from this recursive discrete Fourier operation circuit is 
delayed by a sampling period by the delay unit 42 and then 

35 the signal is multiplied with Ts = sin (2jtk/N) by the multiplier 
34 . In other input signal to be supplied to the aforementioned 
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adder 38, other part of a signal delayed by a sampling period 
by the delay unit 42 is multiplied with Tc = cos(2jtk/N) by 
the multiplier 35. 

Such summed signal by the adder 37 is supplied as real 
5 part signal X r (k, t+1) of an output signal from the recursive 
discrete Fourier operation circuit and part of the signal is 
supplied to the delay circuit 41 as a recursive signal. The 
summed signal by the adder 38 is supplied as imaginary part 
signal X ± (k, t+1) of an output signal from the recursive 

10 discrete Fourier operation circuit and part of this signal 
is supplied to the delay circuit 42 as a recursive signal. 

Here, the multiplication coefficients Tc and Ts given 
to the multipliers 32-35 are set up by the basic frequency 
setting portion 2 shown in Fig.l. In the recursive discrete 

15 Fourier operation circuit , N is set up for an object resolution 
according to this set frequency and the coefficient of the 
discrete Fourier transformation is computed depending on the 
N and then output ted. 

In this way, complex Fourier transformation result X r (k, 

20 t + 1) - JXi(k, t + 1) can be obtained with respect to the supplied 

data stream x(t) , x(t+l) , x(t+2) , x(t+3) x{t+N-l) , x(k, 

t + 1). 

In the above equations (11), (12), with cos(27tk/N) and 
sin(2xtk/N) as Tc and Ts, the recursive DFT operation circuit 
25 shown in Fig. 4 has been described. These Tc and Ts indicate 
the basic frequencies and the recursive DFT operation circuit 
shown in Fig. 4 is a module for carrying out the DFT operation 
to the basic frequency determined by a single k. Thus, by 
connecting N modules in parallel, it is possible to construct 
30 the DFT operation circuit for all N-point basic frequencies. 

Fig. 5 shows the structure of the recursive DFT operation 
device for N points constructed of arithmetic operation modules 
using N basic frequencies. 

Tc and Ts in each module are expressed as follows as 
35 shown in the same Figure. 
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rc-cos^j - (15) 

rj-sin^J - (16) 

The value k in Tc and Ts of each module In the N- point recursive 
DFT operation device is set up to be N different basic 
5 frequencies under Oszk -1 . 

Consequently, the recursive DFT operating portion for 
carrying out the N-point Fourier transformation with parallel 
operations can be achieved. 

Arithmetic operations in the recursive DFT operating 
10 portion achieved in this way are carried out every sampling 
period t according to the equations (11), (12). By changing 
the value k used in that operation from 0 to N-l, a result 
of Fourier transformation under each basic frequency can be 
obtained. 

15 Because the recursive DFT operation executed using these 

transformation equations computes a Fourier transformation 
value at a predetermined specific point , a Fourier 
transformation value can be obtained independently without 
using any Fourier transformation value of other point or mutual 

20 dependence . 

Because the Fourier transformation value at a 
predetermined frequency can be obtained independently, for 
example, if attention is paid to a specific frequency point 
or k is selected to a specific value, the value of an expression 

25 of trigonometric function shown in the aforementioned 
equations (11), (12) is a constant. 

Comparing the conventional DFT operation, because the 
definition of the conventional DFT is carried out by the 
aforementioned equation (1) and butterfly operation in which 

30 the value of trigonometric function indicating the basic 
function changes with the value k is carried out, calculation 
at a specific point is difficult. However, according to the 
recursive DFT operation shown here, because the value expressed 
with the trigonometric function is constant, calculation of 
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a Fourier coefficient for a specific point is easy. 

This recursive DFT can obtain a result of the DFT 
transformation at a selected frequency point with a simple 
calculation by selecting the value k corresponding to a 
5 frequency to be analyzed. 

Further, because the Fourier analysis can be carried 
out by selecting only a basic frequency desired to be subjected 
thereto, the number of points to be handled in a conventional 
FFT is not restricted to points of exponentiation number based 
10 on 2 (value indicated with power of 2; 2, 4, 8, 16, 32, 64, 
128, • • ) and the Fourier analysis can be carried out by 
specifying any point number. 

Next, the structure of the recursive DFT operation 
circuit in which the arithmetic operation modules described 
15 above are disposed in parallel will be described. 

The recursive Fourier transformation operation for N 
points shown in Fig. 5 has components which carry out the same 
operation in the N modules disposed in parallel and thus a 
Fourier transformation circuit including a common portion for 
20 carrying out the operation instead of the plural components 
can be structured. 

Fig. 6 shows the structure of the N-point recursive DFT 
transformation device in which part of the modules is made 
common . 

25 In the DFT shown in the same Figure, of the modules shown 

in Fig. 5, portions carrying out the same operation are made 
common while not common portion are made into sub-modules. 

In the recursive DFT transformation circuit shown in 
Fig . 4 composed of modules , the delay unit 11 which accumulates 

30 data supplied by sampling by the amount corresponding to N 
samples, the subtractor 12, which subtracts the oldest data 
value stored temporarily in the delay unit 11 from the data 
supplied by sampling and the multiplier 31, which multiplies 
a signal supplied from the subtractor 12 with a predetermined 

35 coefficient A, are used as common components. 

The sub-modules 1 to (N-l) are so constructed that such 
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common components are removed from the modules 1 to (N-l). 
Each of the sub-modules 1 to (N-l) utilizes a signal supplied 
from the common delay unit 11, subtractor 12 and multiplier 
31 , so as to reduce the size and power consumption of the circuit . 
5 The recursive discrete Fourier transformation device 

has been described above . Further , as a device , which performs 
transformation of signal format inversely, there is arecursive 
inverse discrete Fourier transformation device. 

Although the recursive discrete Fourier transformation 

10 device is a device for transforming time-series data to 
frequency-basis data, the recursive inverse discrete Fourier 
transformation device is a device for transforming the 
frequency-basis data to time-series data. The recursive 
inverse discrete Fourier transformation device can be realized 

15 by applying the construction method of the device described 
here. 

Then, when the frequency-basis data is transformed to 
time- series data by the recursive inverse discrete Fourier 
transformation device and then, the time- series data is 

20 transformed to the frequency-basis data by the recursive 
discrete Fourier transformation device, if both use 
corresponding parameters for transformation, its original 
frequency-basis data is regenerated. 

Such parameters are numerical values for specifying the 

25 amplitude indicated as A in the aforementioned equations ( 13 ) , 
(14) and as this numerical value, 1, N or square root of N 
is allocated. 

That is , in case where symmetrical signal transformation 
or signal retransf ormation is carried out by supplying data 

30 generated by recursive inverse discrete Fourier 
transformation using, for example, numeric value a as a 
parameter to recursive discrete Fourier transformation, if 
the Fourier transformation is carried out using the parameter 
A corresponding to the parameter a for the Fourier inverse 

35 transformation, data prior to the transformation can be 
obtained. 
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If the amplitude of a signal given with the a and A is 
1 on the Fourier inverse transformation side, 1/N is given 
on the Fourier transformation side, and if 1/N is given on 
the Fourier inverse transformation side, 1 is given on the 
5 Fourier transformation side and if the Fourier inverse 
transformation side depends on an inverse number of square 
root of N, the Fourier transformation side also give an 
amplitude based on an inverse number of square root of N. 
Consequently, data prior to transformation can be obtained 
10 by Fourier-transforming a signal subjected to Fourier inverse 
transformation . 

When this relation has been described with reference 
to the aforementioned equations (11), (12), 
1 

15 indicating the amplitude is assumed to be A and then the 
recursive DFT transformation circuit shown in Fig. 4 has been 
described. This value has been introduced for definition of 
the DFT based on the aforementioned equation and originally 
is a coefficient, which exists for adjustment of the amplitude 

20 between the DFT and the IDFT. 

Therefore, the multiplier 3 is designed such that the 
constant A in the recursive DFT is selected as a value 
corresponding to a constant a in the corresponding IDFT, an 
appropriate value which is in complementary relation with the 

25 IDFT side such as 1 or 1/N as well as 
1 

The Fourier inverse transformation is employed in 
orthogonal frequency division multiplexing modulation (OFDM) 
signal generating unit and the Fourier transformation is 
30 employed in OFDM signal receiving unit so as to generate a 
signal to be digital-modulated and transmitted by the Fourier 
inverse transformation . On the receiving side , a transmitted 
signal is decoded by Fourier transformation and the constant 
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a and A for determining the amplitude may be set up to a 
corresponding value as described above. 

However, because the OFDM signal to be generated and 
transmitted actually in such a manner has many factors which 
5 may change the amplitude value (gain of a signal to be 
transmitted) in its process, the amplitude at the discrete 
Fourier transformation time does not has to accurately coincide 
with a value at the inverse discrete Fourier transformation 
time and thus the value A in discrete Fourier transformation 

10 may be set up by selecting a value easy to be incorporated 
appropriately . 

The discrete Fourier transformation introduced from the 
equations (11), (12) and the inverse discrete Fourier 
transformation as its application have been described. 

15 These equations for obtaining the Fourier coefficient 

is capable of obtaining a Fourier coefficient value in a simple 
form by using a difference between two data values separate 
by N samples, and sine and cosine of a basic frequency with 
respect to a Fourier coefficient obtained before a sample 

20 clock. 

The cosine Tc and sine Ts of the basic frequency used 
in the discrete Fourier transformation introduced from the 
equations (13), (14) indicate that the recursive DFT 
transformation device can be constructed by using functions 

25 orthogonal to each other, and thus by using orthogonal 
functions to each other, for example, replacing the Tc with 
sine Ts, the same analysis method can be realized. 

For the purpose, a combination of a relation between 
the real part and imaginary part of a Fourier coefficient for 

30 use as a recursive signal and the polarities of these signals 
needs to be determined so that the Fourier coefficient can 
be obtained in as a simple format as possible and then , a Fourier 
coefficient has to be introduced according to that determined 
method . 

35 < Second Embodiment > 

Next, a second embodiment of the recursive discrete 
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Fourier transformation device of the present invention will 
be described. In a following description, a description of 
common components as the first embodiment is omitted. 

Fig . 7 shows a second embodiment of the recursive discrete 
5 Fourier transformation device and it is described with the 
Figure. 

The recursive discrete Fourier transformation device 
comprises a data updating portion 101 to which data sampled 
in a predetermined time interval is supplied and which stores 

10 temporarily supplied new N (N: positive integer ) data, a basic 
frequency setting portion 102 in which the basic frequency 
for carrying out the discrete Fourier transformation is set 
up, a recursive processing portion 103 for carrying out signal 
processing while using a signal temporarily stored recursively 

15 and a multiplication processing portion 104 for carrying out 
multiplication processing and outputting an obtained Fourier 
transformation coefficient. 

Next, an operation of the recursive discrete Fourier 
transformation device having such a structure will be 

20 described. 

An operation of the data updating portion 101 is the 
same as the first embodiment. 

Next, N data stored temporarily in the data updating 
portion 101 is supplied to the recursive processing portion 

25 103 and thus Fourier transformation operation is carried out 
by the recursive processing portion 103 and the multiplication 
processing portion 104 according to a method, which will be 
described later. 

The recursive processing portion 103 stores an 

30 intermediate result of arithmetic operation temporarily in 
a memory circuit contained therein and carries out signal 
pre-processing for executing Fourier arithmetic operation 
while using the temporarily stored signal recursively. The 
pre-processed signal is supplied to the multiplication 

35 processing portion 104 , in which the multiplication processing 
as post processing is carried out so as to obtain the real 
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part X r (k, t + 1) and imaginary part Xi(k, t+1) of the Fourier 
coefficient . 

In order to obtain a result of Fourier arithmetic 
operation analyzed with a predetermined resolution, a basic 
5 frequency for giving a desired resolution is set up by the 
basic frequency setting portion 102 and operation parameters 
of the recursive processing portion 103 and the multiplication 
processing portion 104 are set up. Following the set 
parameters, the result of the Fourier analysis based on 

10 predetermined frequency resolution information is computed 
and outputted. 

Next, this recursive discrete Fourier operation method 
will be described in detail with reference to the Fourier 
operation result conventionally carried out . 

15 In this second embodiment , the equations (11), (12) are 

introduced like the first embodiment so that complex Fourier 
transformation result X r (k, t+l)-jXi(k, t+1) can be obtained 
at time t+1. 

It means that a result of Fourier transformation on the 
20 newest N-point data containing new data obtained by successive 
sampling at the sampling time interval shown in Fig. 2 can be 
obtained. 

Here, if it is regarded that the Fourier transformation 
executed in the above-describedmanner is a linearly invariable 

25 system in which continuous sampling data is supplied to a system 
having a certain coefficient and converted, this conversion 
processing can be regarded to be filtering processing. 

Then, if l/( square root of N) which gives an amplitude 
is assumed to be A and the real part X r (k,t) and the imaginary 

30 part X± ( k , t ) of Fourier coefficient with respect to data stream 
supplied since time t are assumed to be Y r , k (t) and Y i<k (t), 
the aforementioned equations (11), (12) are treated as a 
transformation equation of output y to input x, so that the 
transformation equation of the Fourier coefficient of data 

35 supplied since time t+1 can be expressed as follows. 
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3^(f+l)-{y,*W + 4*(' + tf)-^^^ - (17) 

yJt + l)-y„(t)cos[2fyb,t(t) + 4^ .» (18) 



( where A = —== ) 

4n 



5 Further, if these equations (17), (18) are subjected 

to Z- transformation, a result of the transformation can be 
expressed as shown in following equations (19), (20). 

n, M = fc, (*V l + 4x(zV - *{*)}* cos [ 2 ^] + M* -1 sin [ 2: f ] - ( 19 ) 

raM-ra^co^^ (20) 

10 If a transfer function of output Y(z) to input X(z) is 

assumed to be H(z), H(z) is given in a following equation. 

Then, if the equations (19), (20) are treated as 
simultaneous equations and it is intended to obtain respective 
15 transfer functions of the real part and imaginary part, they 
can be obtained as shown in following equations (21), (22). 



hXz) = a{i-z-») 



\l-2cos\2^\z' 1 +z- : 



H i {z) = A(l-z- N ] 



l-2cos 2 — \z~ l +z~ 
N\ 



- (21) 



(22) 
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Like this, a transfer function indicating the recursive 
DFT is expressed in a following equation (23). 



H(z) = a(1-z- n 



cos [2— 1 - j sin [2— 1 

1 *l I *M 



l-2cos 2 — z" 1 +z" : 
N\ 



(23) 



In this equation, 1-Z" N indicates deletion of the oldest 
5 data upon updating the input data and a portion expressed with 
a fraction is so constructed that its numerator and its 
denominator are IIR filters (cyclic filter) having primary 
and secondary degrees respectively. 

Although the DFT is a FIR filter for carrying out 
10 processing to a finite point number N. if it is achieved with 
a IIR filter as shown in the aforementioned equation (23), 
it can be constructed with a simpler hardware than the FIR 
filter. 

Next, a method for achieving the characteristic 
15 expressed with this transfer function with hardware will be 
described. 

Here, if the basic frequencies cos (2jik/N) and sin( 2jtk/N) 
are assumed to be Tc, Ts (gamma c, gamma s) , the equation (23) 
can be expressed as follows. 

20 *W-4-^i^=Srl 



(Where. Tc - cos [ 2 f ] . O-smJzf]. A-j=, 

A method for structuring the recursive DFT based on the 
equation (24) will be described. 

Fig. 8 shows a recursive DFT operation unit for carrying 
out the discrete Fourier transformation and the structure and 
operation thereof will be described in detail. 

The recursive DFT operation unit comprises an updating 
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processing portion 101 including a memory 111 and a subtracter 
112, a recursive processing portion 103 including a subtractor 
131, an adder 132, memories 133, 134 and a multiplier 135 and 
a multiplication processing portion 104 including multipliers 
5 141, 142, 143 and 144. 

In the same Figure, the memory 111 is notated with Z" N 
and the memories 133, 134 are notated with Z" 1 . The number 
of this exponent indicates a period during which supplied data 
is stored temporarily. -N indicates temporarily storing in 

10 N sampling period while -1 indicates temporarily storing in 
a sampling period. 

Further, the subtractors 112, 131 subtracts signals 
supplied in a vertical direction ( from down to up ) from signals 
supplied in a horizontal direction in the same Figure and a 

15 signal obtained by the subtraction is supplied as subtraction 
output . 

The adders 132 , 144 sum up two supplied signals indicated 
with arrows, so that a signal obtained by the summing is 
outputted . 

20 The multipliers 135, 141, 142, 143 multiply a supplied 

signal with a value based on a symbol indicated on the side 
of a multiplier in the same Figure and supply a result of the 
multiplication. The values with which those multipliers 
multiply are -2rc, Arc, -A, -Ars. 

25 Next , an operation of the recursive DFT operation device 

having such a structure will be described. 

First, x(t+N) is supplied to the subtractor 112 and the 
memory 111 as an input signal. 

The memory 111 delays a supplied signal by N sampling 

30 periods. The output signal of the memory 111 at this time 
is supplied as a signal supplied before the N sampling periods , 
that is , x( t ) is outputted while it is supplied as a subtraction 
signal to the subtractor 112. Thus, the updating processing 
portion 101 carries out arithmetic operation of x(t+N)-x(t) . 

35 The updating processing portion 101 updates data 

accumulated in each of N sample clocks . In another view, such 
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signal processing is a processing in which an input signal 
is multiplied with 1-Z" N which is a first transfer function. 

A signal obtained by such processing is supplied to the 
subtractor 131 in the recursive processing portion 103. Here, 
5 a signal from the adder 132 is subtracted and part of signals 
obtained by the subtraction is supplied to the multiplication 
processing portion 104 while the other part thereof is supplied 
to the memory 133. 

Signals delayed by a sample period by the memory 133 

10 are supplied to the multiplier 135, the memory 134 and the 
multiplier 142 of the multiplication processing portion 104. 

A signal multiplied with a predetermined number by the 
multiplier 135, supplied by the memory 133 is supplied to an 
input terminal of the adder 132. Then, a signal delayed by 

15 one more sample period by the memory 134 is supplied to other 
input terminal of the adder 132 and summed up with the 
aforementioned signal, and the summed signal is supplied to 
a subtraction signal terminal of the subtractor 131. 

In other words, a signal obtained by providing a signal 

20 outputted from the subtractor 131 with the characteristic - 2T C 
XZ" 1 from the memory 133 and the multiplier 135 and a signal 
obtained by providing a signal outputted from the subtractor 
131 with the characteristic Z~ 2 from the. memories 133, 134 
are summed up by the adder 132 . Then , the result of the summing 

25 is supplied to the subtraction signal terminal of the 
subtractor 131. The subtractor 131 subtracts a result of the 
summing from a signal from the updating processing portion 
101. Therefore, the recursive processing portion 103 carries 
out multiplication processing expressed with a second transfer 

30 function l/(l-2r c Xz" 1 + Z" 2 ) to carry out arithmetic operation 
using old data recursively. 

A signal ql obtained by multiplying with a second 
transfer function and a signal q2 obtained by delaying the 
signal ql by a sample period by the memory 133 are supplied 

35 to the multiplying portion 104. 

In the multiplication processing portion 104, one of 
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the signals ql is multiplied with Arc by the multiplier 41 
to obtain a first multiplication signal and the signal q2 is 
multiplied with -A by the multiplier 142 to obtain a second 
multiplication signal. These first and second multiplication 
5 signals are supplied to the adder 144 , in which they are summed 
up to obtain a signal X r (k, t+1). This obtained signal is 
supplied as a real-part signal of a Fourier coefficient. 

The other of the signals ql supplied to the 
multiplication processing portion 104 is multiplied with -Arc 

10 to obtain a signal Xi(k, t+1) and the obtained signal is supplied 
as an imaginary part signal of the Fourier coefficient. 

In arithmetic operation of the transfer function in the 
multiplication processing portion 104, an operation for the 
real part and an operation for the imaginary part with respect 

15 to a third transfer function A(T C - jr s -Z _1 ) are executed 
independent of each other so as to obtain a complex Fourier 
coefficient output signal. 

As described above, the recursive DFT arithmetic 
operation device can introduce complex Fourier transformation 

20 result X r (k, t + l)-jX i (k / t+1) at time t+1 by using data x(t+N) 
sampled and supplied at time t+1 and using recursively 
intermediate data under Fourier operation for data supplied 
before . 

Although in Fig. 8, the signals supplied from the 
25 recursive processing portion 103 are two signals, ql and q2, 
if the multiplication processing portion 104 contains the 
memory 133, it is permissible to omit a signal line of the 
q2. 

Fig. 9 shows the structure in case where the recursive 
30 processing portion 103 and the multiplication processing 
portion 104 are connected through a single signal line ql. 

However, the multiplication processing portion 104 
requires additionally a memory 145 compared with the 
multiplication processing 104 of Fig. 8. 
35 In the above descriptions, the recursive DFT operation 

device can be constructed by dividing the transfer function 
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shown in the aforementioned equation (24) to first transfer 
function 1-Z" N , second transfer function 1/ ( l-2r c Z" 1 +Z" 2 ) and 
third transfer function A(r c -jr s -Z" 1 ) . By constructing the 
arithmetic logic circuits having these delay, 
5 addition/subtraction and multiplication functions with 
digital circuits, it is possible to configure the recursive 
discrete Fourier transformation circuit with a smaller, simple 
circuit block. 

According to the recursive discrete Fourier 

10 transformation method using the transfer function, the Fourier 
analysis can be carried out by selecting a desired basic 
frequency. Thus, the number of points to be handled is not 
restricted to power of 2 unlike in case of the FFT, so that 
any positive integer can be adopted. 

15 Further , in the Fourier transformation carried out using 

any integer N, it is possible to perform a frequency analysis 
with high resolution in time direction by changing the length 
N of supplied data stream, for example, increasing the value 
of N . That is , because an operation result , in which the degree 

20 of analysis resolution of Fourier transformation is set to 
a desired value, can be obtained by changing the value N, the 
recursive discrete Fourier transformation circuit having such 
a structure can execute the Fourier transformation processing 
by selecting any resolution. 

25 As described above, the recursive DFT operation device 

can execute the Fourier transformation processing by selecting 
any resolution. The frequency resolution can be set to any 
resolution depending on each frequency subjected to analysis . 
According to the second embodiment, the Fourier 

30 transformation can be carried out to individual frequencies 
at the same time with the structure shown in Fig. 5 like the 
first embodiment. 

Further, it is permissible to construct a N-point 
recursive DFT transformation device whose part of the modules 

35 are used in common like the first embodiment. 

Fig. 10 shows the structure of the N-point recursive DFT 
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transformation device whose part of the modules are used in 
common according to the second embodiment. 

In the recursive DFT transformation circuit shown in 
Fig. 8, composed of modules, the delay units 111 for 
5 accumulating data sampled and supplied by the amount 
corresponding to N samples and the subtracters 112 which 
subtract the oldest data stored temporarily in the delay unit 
111 from data value sampled and supplied are used in common. 
The sub -modules 1 to (N-l) are so constructed that these 

10 common components are removed from the modules 1 to (N-l). 
The sub-modules 1 to (N-l) employ signals supplied from the 
common delay units 111 and subtracters 112 in common, thereby 
achieving reductions of circuit size and power consumption. 
According to the second embodiment, the recursive 

15 inverse discrete Fourier transformation device , which carries 
out transformation of signal format inversely, can be 
constructed like the first embodiment. Therefore, a 
description about the recursive inverse discrete Fourier 
transformation device of the first embodiment is applicable 

20 to the second embodiment just as it is. 

It should be understood that many modifications and 
adaptations of the invention will become apparent to those 
skilled in the art and it is intended to encompass such obvious 
modifications and changes in the scope of the claims appended 

25 hereto. 
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